Interpretation of low-temperature nuclear quadrupole resonance spectra in 
Lai 875Bao.i25Cu04 in terms of two-dimensional spin superstructure 
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This paper reanalyzes the low temperature nuclear quadrupole resonance (NQR) experiments in 
Lai.875Bao.i25Cu04 by Hunt et al. [Phys. Rev. B 64, 134565 (2001)] in an attempt to determine 
the dimensionality of spin modulations in this and other compounds of the lanthanum family of 
high temperature cuprate superconductors. It is concluded that the shape of the NQR spectra 
obtained by Hunt et al. favors the two-dimensional pattern of spin modulations known as "grid" or 
"checkerboard". The paper also contains the discussion of charge patterns, which can accompany 
the above spin modulation. 
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I. INTRODUCTION 

The notion of stripes has become a standard concept 
in the field of high-temperature superconductivity about 
ten years ago, after Tranquada et al.— have observed the 
four-fold splitting of the antiferromagnetic (AF) (tt, tt) 
peak in elastic neutron scattering experiments in the lan- 
thanum family of cuprates. The above authors have in- 
terpreted the resulting four peaks, as coming from two 
different kinds of stripe domains, each exhibiting a unidi- 
rectional spin modulation (referred to below as ID) run- 
ning along one of the two principal lattice directions. The 
ID spin modulation was to be accompanied by the ID 
charge modulation peaked around the antiphase mag- 
netic boundaries, hence the name "stripes" . There ex- 
ists, however, an alternative interpretation sometimes re- 
ferred to as grid or checkerboard, which stipulates that 
the four-fold splitting of the (tt, tt) neutron peak is caused 
by a two-dimensional (2D) modulation of the AF order 
shown in Fig. [T] Such a modulation is accompanied by 
a 2D charge density wave localized around the antiphase 
boundaries, which have the appearance of grid — the 
term to be used in the rest of this paper. 

Grid superstructure has been mentioned in the litera- 
ture since the early numerical works on the subjec1>2ii2ii^i^, 
however, it was only quite recentlj*^ii, that grid was ad- 
vocated as a preferred interpretation of the (tt, tt) peak 
splitting. It has also been suggested in Ref.^ that similar 
2D inhomogeneous structures are universally present in 
other families of cuprates in the normal and supercon- 
ducting states. The insight of that work was not com- 
pletely trivial, since its appearance predated the scanning 
tunneling spectroscopy (STS) studies of Vershinin et al^, 
Hanaguri et al^, and the neutron study of Hinkov et ali^ 
— all three are the works, that brought the dimensional- 
ity of generic spin modulations in cuprates back into the 
spotlight (see, e.g.fi^). 

The connection between the inhomogeneous structures 
in lanthanum cuprates and the generic properties of other 
high-temperature cuprate superconductors is a delicate 
issue. Spin and charge modulations observed in this 




FIG. 1: [Color online] Commensurate diagonal grid modula- 
tion [cite-centered] . The antiferromagnetic order is modulated 
by function ((2|. The length of the arrows is proportional to 
the value of staggered magnetization. Lines drawn through 
the zero magnetization sites indicate antiphase boundaries of 
magnetic domains. They form a "grid" pattern. 



cuprate family may turn out to be a composition-specific 
feature, or, in the opposite extreme, similar modulations 
may be universally present in all families of cuprates 
but observable only in the lanthanum family due to 
anomalous slowing down. In this context, the differ- 
ence in the dimensionality of the modulation is not in- 
consequential. The apparent message contained in the 
symmetry of the superconducting order parameter—, the 
symmetry of the pseudogap^'^, the STS observations on 
the surface of Bi2Sr2CaCu208+5 (Bi-2212)^J^ii^ and 
Ca2_a;Na2;Cu02Cl2^ is two-dimensional. Therefore, if 
the ID nature of the spin modulations in lanthanum 
cuprates is established conclusively, then, in the opinion 
of this author, it will likely imply the composition-specific 
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nature of the phenomenon. (The authors of Refsj^^ii^ 
seem to have the opposite expectation.) In contrast, 
if the spin modulations in this cuprate family are two- 
dimensional, then this would be yet another indication 
that the tendency to form 2D modulations is a general 
property of cuprates and, therefore, may be essential for 
the mechanism of superconductivity. The differences be- 
tween different geometries of 2D modulations should be 
much less important. In lanthanum cuprates, magnetic 
modulations may develop into a diagonal grid pattern, 
while, in other cuprates, the system may phase separate 
into more-AF-correlated and less-AF-correlated nano- 
clustersi^ without forming quasi-static magnetic order 
and, therefore, not accompanied by antiphase bound- 
aries. 

Since the argument for the non-universality of the 
modulations in lanthanum cuprates has merit, it is im- 
portant to establish the dimensionality of these modula- 
tions without invoking the knowledge about other fam- 
ilies of cuprates or nickelates. The discussion in the 
literature specific to lanthanum cuprates has been lim- 
ited so far mostly to ncutron^-iii^i2ii2^i2^ and hard X- 
ray^l experiments. The primary focus of this paper is 
on extracting the dimensionality of spin modulations in 
Lai.875Bao.i25Cu04 from the nuclear quadrupole reso- 
nance (NQR) results of Hunt et al^ in combination with 
the muon spin rotation (/xSR) results of Nachumi et alM-. 
The paper also elaborates on some of the dimensionality- 
related arguments made earlier in Ref.— . Most of the pa- 
per will contrast ID stripes with 2D grid, but ID spin 
spirals22i22iiSS will also be discussed in Section |Vl 

The plan of the rest of the paper is as follows. In 
Section [Til the present non-NQR knowledge about the 
dimensionality of spin and charge modulations in lan- 
thanum cuprates is discussed. Section IIIII discusses rel- 
evant NQR findings of Hunt et al. The NQR-related 
findings of the present work are presented in Section HVl 
and discussed in Section fVl 

Sections HIl and [TTTl amount [mainly] to an extended in- 
troductiry review. The readers familiar with the broader 
issues involved can skip the above sections and proceed 
directly to Section HVl 



II. CHARGE MODULATIONS IN THE 
LANTHANUM FAMILY OF CUPRATES 

Magnetic peaks corresponding to higher order modu- 
lation harmonics have not been observed so far and thus 
cannot be used to discriminate between different mod- 
ulation patterns. A direct resonant X-ray observation 
of charge carrier density modulations reported recently 
by Abbamonte et al.'^'^ has not addressed the dimen- 
sionality issue either. The direct attempts to discrimi- 
nate between ID stripe and 2D grid patterns have been 
made so far only on the basis of the orientation of charge 
satellites^iiiii^, which represent the response of the lat- 
tice to the periodic charge and/or spin modulations. 



In Refsiiiii^, the analysis of experiments is based on 
Landau expansion^! , which requires small spin and 
charge modulations. The resulting charge patterns for 
the ID and the 2D cases are shown in Fig. [51Ja,b). In the 
framework of Landau expansion, leading charge peaks 
corresponding to ID stripe and 2D grid pictures do not 
coincide, and, therefore, the observation of charge satel- 
lites expected for the ID case together with the non- 
observation of the satellites expected for the 2D case 
(notwithstanding the weakness of the charge signal) do 
present a clear argument. 

There are, however, good reasons to expect that the 
amplitude of the spin modulations is not small. For ex- 
ample, as indicated by /zSR^^, the amplitude of the mod- 
ulated staggered magnetization in Lai.875Bao.i25Cu04 is 
large — about half of the staggered moment of the un- 
doped parent compound. Strong modulation of staggered 
magnetization would induce a deeper local potential for 
charge carriers thus forcing them to be stronger localized 
around the antiphase boundaries. Larger amplitude of 
charge modulation implies stronger interaction between 
the charge carriers, which would undermine the theoret- 
ical basis for the Landau expansion and introduce new 
charge peaks. One particular problem with the Landau 
expansion is that it gives the maximum of charge density 
around the nodes of the antiphase grid [see Fig. (Hb)]. 
If the charge modulations are large, then the Coulomb 
repulsion between the charge carriers would rather lead 
to a smaller charge density around the nodes. 

In the case of large charge modulations, the leading 
charge peaks corresponding to the ID and 2D interpre- 
tations of the magnetic (tt, tt) peak splitting can coincide. 
If charge density is reasonably localized around the an- 
tiphase grid but avoids grid nodes, then a useful insight 
can be gained from analysing the density profile shown 
in Fig.[2l^c). It is given by function 
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where x and y are the coordinates along two principal 
lattice directions and a is the lattice period. This profile 
was mentioned in Refi^ but not elaborated upon. It cor- 
responds to the same splitting of magnetic (f , f ) peak 
as the one arising from the ID stripes spaced by 4a. The 
cos^ factors in Eq.([T|) set the grid pattern, while the 
sin" factors suppress the intensity near the grid nodes. 
The positions of the resulting charge peaks and their 
relative intensities are shown in Fig. [3l (See the Ap- 
pendix for the calculation.) Leading charge peaks cor- 
responding to ID stripes spaced by 4a and to modu- 
lation (III) coincide with each other and are located at 
(±^,0) and (0, ±^), where charge peaks were detected 
in Refsiii^i^i22ii2i. Modulation ([T]) also generates peaks 
at (±^,±^) predicted by Landau expansion, but these 
peaks have only third largest intensity (about 7 times 
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FIG. 2: [Color online] Three color-coded patterns of charge 
modulations corresponding to the same splitting of magnetic 
(tt, tt) peak: (a) stripes; (b) grid in the Landau approxima- 
tion; (c) grid with the charge density avoiding the grid nodes 
[Eq.([l|]. Dashed lines represent antiphase boundaries of an- 
tiferromagnetic domains. 



FIG. 3: [Color online] Positions of charge peaks (filled circles) 
corresponding to the scattering from the charge modulation 
given by Eq.(IT| and shown in Fig. O^c). Numbers above the 
circles represent the relative intensity of the peaks obtained 
in the Appendix. Large empty circles and empty diamonds 
are the positions of leading charge peaks corresponding, re- 
spectively, to stripes [Fig. [IJa)] and to grid in the Landau 
approximation [Fig. [D^b)] 



smaller than the intensities of the leading peaks). The 
peaks of the second largest intensity (4 times smaller than 
that of the leading peaks) are located at (±^,±^). 

It may appear that the leading peak position for the 
density profile llj is the result of a fine tuning of the 
powers of sine and cosine. This is, however, not the case. 
The only important aspect of this profile is that the fac- 
tors sin^ suppress the charge density near the grid nodes, 
as required by the Coulomb repulsion between charges. 
As illustrated in the Appendix, if the sin^ factors are re- 
moved, then the leading peaks are located at the position 
predicted by the Landau expansion, but, if they are left, 
then the position of the leading peaks at (±^,0) and 
(0, ±2^) does not depend on the power of the cosine fac- 
tors. The unrealistic aspect of profile H]) is that the sin^ 
factors suppress the charge density near the grid nodes 
completely. However, even an incomplete Coulomb sup- 
pression of the charge density near the grid nodes can 
easily shift the intensity maximum from the Landau ex- 
pansion result to (±^,0) and (0,±^). 

From the above perspective, it is the non-observation 
of the weaker (±^,±4^) and (±t?-,±t|^) peaks rather 
than the observation of the (0,±^) peaks that would 
be the real argument in favor of the ID interpretation. 
However, to the best of this author's knowledge, the open 
literature contains only one experimental scan across the 
ia) ^'^^ (ia' ia)' "^^ich has been presented in Refi^ 
for Lai.48Ndo.4Sro.i2Cu04. In that case, the experimen- 
tal resolution was just not sufficient to claim observation 
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or non-observation of the peaks, which would be 4 or 7 
times weaker than the leading peak observed at (0, ^). 
In fact, quite remarkably the scan around (5^, 5^) shows 
a promising sign of a peak. Repeating that scan with a 
better resolution would certainly have merit. 

Another argument frequently cited in support of the 
ID stripe picture is that the observability of static spin 
and charge modulations is correlated with the onset 
of the "non-generic" low temperature tetragonal (LTT) 
crystal phase. Tranquada et al^ have proposed, that 
the ID stripes follow the lattice anisotropy of the LTT 
phase. This anisotropy alternates between adjacent 
Cu02 planes, which would imply, that stripes in two 
adjacent planes run perpendicular to each other. Em- 
pirically, however, the onset of the LTT phase is not a 
necessary condition for the observability of static mod- 
ulations, because quasielastic signatures of spin modula- 
tions have also been reported for the "generic" low tem- 
perature orthorhombic (LTO) phase of La2-xSr2:Cu04 
with and without Zn codoping22ii2^i^. But one can still 
hope that the presence of the LTT phase is at least a 
sufficient condition for the observability of static modu- 
lations. 

Although the above proposition of Tranquada et al. 
has visual appeal, the onset of the LTT phase as such 
does not readily discriminate between the stripe and the 
grid interpretations. The general expectation is that the 
magnetic (tt, tt) peak splittings observed by elastic neu- 
tron scattering in the LTT phase and by inelastic neutron 
scattering in the LTO phase have the same origin. If true, 
this implies that the system has chosen the dimension- 
ality of the modulation pattern (ID or 2D) without the 
LTT phase being involved, and then the LTT phase only 
plays a role in slowing the modulations down, to make 
them observable in elastic experiments^Si. Therefore, to 
claim that the symmetry of the LTT phase discriminates 
between the grid and the stripe pictures is equivalent to 
claiming that the LTT phase would slow down fluctu- 
ating stripes, but would not slow down fluctuating grid. 
This author is not aware of any experimental or theoret- 
ical result, which would support the latter part of this 
claim. 

The picture of stripes running in perpendicular direc- 
tions in the adjacent planes has further difficulty that 
such a configuration is not most favorable energetically 
from the viewpoint of interplane Coulomb interaction. 
The lowest Coulomb energy corresponds to stripes run- 
ning in the same direction in two adjacent planes with 
pattern in one plane shifted with respect to the other by 
half a period of charge modulation. Since the stripes are 
unscreened lines of charge, the interplane Coulomb en- 
ergy should be large. Can the coupling of stripes with 
the lattice anisotropy outweigh the Coulomb coupling be- 
tween stripes from adjacent planes? It does not seem very 
likely. 

Recently, Robertson et aL™ have also argued that, due 
to the anisotropy of the LTT phase, the incommensurate 
grid modulation should have slightly different periodicity 



in the adjacent Cu02 planes, which should lead to 8-fold 
splitting of the (tt, tt) peak instead of the 4-fold splitting 
observed experimentally. Since the crystal lattice peri- 
ods cannot be detectably different in the adjacent planes 
(because of the prohibitive elastic energy cost), the ar- 
gument of Roberson et al. is not very direct, as it relies 
on the additional theoretical assumption, that the mech- 
anism similar to the Fermi-surface nesting is behind the 
formation of the grid superstructure. This argument has 
further difficulty, that the mismatch of the charge super- 
structure periods would lead to the regions, where the 
maxima of charge modulation density in two adjacent 
planes lie on the top of each other, which, in turn, leads 
to a large Coulomb energy cost. 

The critical discussion so far does not change the fact 
that the observed positions of the charge peaks, and 
the involvement of the LTO-LTT transition cannot be 
turned into an argument favoring the grid over the stripe 
pictures, and to this extent, these two observations do 
support the ID stripe interpretation. In the present 
case, however, the the presumption of innocence can 
be rephrased as follows: A two-dimensional experimen- 
tal peak pattern should be considered as an evidence 
for two-dimensional superstructure unless conclusively 
proven otherwise. In the opinion of this author, the above 
cited facts fall far short of delivering a conclusive proof. 

One can, however, start from the "presumption of 
guilt" by considering the 2D grid interpretation of the 
2D neutron data as "exotic", and requiring compelling 
reasons to choose it over the ID stripe picture. It is the 
purpose of the present paper to analyse the NQR ex- 
periments from this perspective, and to point out that, 
within the limits of experimental uncertainties, they do 
constitute a clear indication in favor of the 2D grid in- 
terpretation. 



III. DISCUSSION OF NQR RESULTS BY HUNT 
ET AL. 

A slowly fiuctuating spin modulation, ID or 2D, is 
expected to freeze at sufficiently low temperatures thus 
leading to a broad NQR spectrum having shape, which 
may discriminate between different dimensionalities of 
the modulation. In Ref.[l], Hunt et al. have presented, 
among other findings, the results of their NQR study of 
Lai.875Bao.i25Cu04 aimed at testing the stripe interpre- 
tation of neutron scattering experiments. At tempera- 
tures below 8K, they have, indeed observed a broad NQR 
spectrum of ®^Cu [shown in Fig.Uff)], which signified the 
expected freezing of local hyperfine fields affecting cop- 
per nuclear spins. This study was done in the context of 
a broader effort to elucidate the dynamics of electronic 
spins in lanthanum cuprates by the means of nuclear 
magnetic and quadrupole resonance o'^^'^^i'^^'^^i'^^'^^i"^^ 

The above NQR spectrum has exhibited one clear fea- 
ture, the 31 MHz peak, but an important theoretical find- 
ing of Hunt et al. was that this particular peak should 
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FIG. 4: [Reproduced from Ref.— ] Distributions of hyperfine 
fields and corresponding NQR spectra obtained in Ref.—. 
Solid lines represent theoretical calculations. Circles in fig- 
ure (f) represent experimental NQR spectrum at 350 mK. 
Vertical axes have arbitrary units everywhere. 



be present for any continuous distribution of hyperfine 
fields passing the interval 0.5-1.5 T, where the minima 
of two NQR frequency branches are located - both cor- 
responding to frequency 31 MHz (see Fig. O. Therefore, 
the presence of the 31 MHz peak could not discriminate 
between different kinds of modulations. Apart from this 
peak, the spectrum has exhibited no other structure. 

Hunt et al. have attempted to fit the NQR spectrum 
to the predictions of many possible ID static stripe pic- 
tures, but have found that all such pictures lead to very 
structured lineshapes having several peaks not observed 
experimentally [see Figs. HJa-d)]. 

The origin of the structure in theoretical NQR line- 
shapes is two-fold, (i) For commensurate superstruc- 
tures, there is only a finite number of distinct nuclear 
sites. The traces of this discreetness then appear in the 
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FIG. 5: [Reproduced from Ref.^] "^^Cu NQR transition fre- 
quencies as a function of hyperfine field. The transition fre- 
quencies are obtained from Hamiltonian ((SJ. Dashed lines 
indicate the regions of low probability transitions. 



NQR lineshape. (ii) For uniform incommensurate ID 
modulations, which contain the infinite number of dis- 
tinct sites, the peaks in the NQR lineshapes result from a 
Van Hove singularity in the distribution of ID-modulated 
hyperfine fields. This singularity is located at the maxi- 
mum in the absolute value of the hyperfine fields [see the 
middle panel of FigUI^d)]. 

A priori, a 2D grid modulation represents a promising 
alternative to interpret the rather featureless NQR spec- 
trum of Hunt et al, because, when the grid is commensu- 
rate, it contains a larger number of inequivalent nuclear 
sites, and thus tends to produce more smeared structure 
in NQR lineshape, and, when the grid is incommensu- 
rate, it would not lead to a Van Hove singularity at the 
maximum value of 2D modulated hyperfine fields. (It 
would lead to a singularity around the zero value of hy- 
perfine field, but that singularity is less important. See 
the next section.) 

Hunt et al. have considered the possibility of grid mod- 
ulation, but only the commensurate version with the ori- 
entation rotated by 45 degrees in comparison to the one 
expected from the positions of neutron peaks. In this 
way, they have, nevertheless, obtained a much better 
agreement with experiment in comparison with any of the 
ID interpretations [see Fig. [D^e)], which is not very sur- 
prising in view of the preceeding discussion. Yet, the grid 
pattern used by Hunt et al. generated too much struc- 
ture, which led them to abandon the grid interpretation 
all together and to test another one in terms of the ran- 
dom distribution of hypefine fields. The latter interpreta- 
tion gave nearly perfect agreement with the experimental 
NQR spectrum [Fig. [Iff)]. The problem with that inter- 
pretation is that the random configuration of staggered 
spin polarizations is inconsistent with the modulations 
of AF spin structure observed by neutron scattering^. 
Random picture should also lead to a monotonic decay 
of the fiSH asymmetry signal rather than the oscillations 
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observed in experiment-^. 

The frequency of the /j,SR asymmetry oscillations im- 
plies that the modulation amplitude of the staggered po- 
larization is about where /is is Bohr magneton. 
This number imposes a difficult constraint on the inter- 
pretations of NQR spectrum. Hunt et al, have concluded 
that any kind of static broadening of the hyperfine field 
distribution corresponding to the modulation amplitude 
0.3/Lt_B would fail to reproduce their NQR spectrum, be- 
cause it would introduce too much spectral weight in the 
higher frequency interval, where much smaller weight is 
observed experimentally. They have found, however, that 
with the modulation amplitude 0.1 — 0.2^b and a very 
significant degree of broadening (50-100 per cent of the 
modulation amplitude), the observed NQR spectra can 
be reproduced. 

Hunt et al. have suggested that the above discrepancy 
may be due to spin fluctuations on the NQR timescale of 
20 fis, which are not seen on the faster ^SR timescale of 
0.5^s. This suggestion, however, runs into the difficulty, 
that, in order to change the shape of the effective distri- 
bution of hyperfine fields, spins should fluctuate on the 
time scale of inverse NQR frequency (~ 0.03 /is) or faster. 
(This, presumably, happens at temperatures higher than 
30 K.) The fluctuations on the scale of 20 /is can lead 
to the loss of resonant spins, which, in turn, can reduce 
the integrated signal intensity but should not change the 
shape of the spectral line. This is consistent with the 
fact, that the NQR spectra of Hunt et al. have identi- 
cal shape but different intensity at 1.7 K and 350 mK. 
Even if, one insists that the spin fluctuations may affect 
the observed NQR lineshape despite the above argument, 
one would have to make a further assumption that these 
fluctuations remain temperature- independent as the tem- 
perature is reduced by a factor of 5 from 1.7 K to 350 K, 
which is quite unlikely. 

The present work has been motivated by the observa- 
tion that Hunt et al. did not consider incommensurate 
grid modulation, which has promise of describing their 
experiment without difficult assumptions. 



IV. NQR SPECTRA FOR DIAGONAL GRID 
MODULATION 

As mentioned in the previous section. Hunt et al. have 
considered grid modulation, which had different orien- 
tation, in comparison to the one that follows from the 
neutron experiments. Their grid was oriented along the 
principal lattice directions, whereas the geometry of the 
neutron (tt, tt) peak splitting implies the diagonal grid 
orientation^ given by the modulation function [leading 
harmonic! 
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and shown in Fig[TJ Here i5 is the dimensionless measure 
of the splitting of the (tt, tt) peak. The four peaks are 



located at [^(1 ± 2(5), f ] and [f , f (1 ± 25)]. The exper- 
iments indicate^ that 5 « i, which corresponds to the 
distance 4a-\/2 between two nearest grid nodes. 

In this section, three cases are computed theoretically 
and compared with experiment: (a) commensurate grid 
corresponding to 5 = |; (b) uniform incommensurate 
grid corresponding to 5 = | — e, where e is an irrational 
number much smaller than i; and (c) disordered incom- 
mensurate grid with 5 « i. Below, unless explicitly in- 
dicated otherwise, the same computational routine and 
the same values of parameters as those of Hunt et al^ 
have been used. 

When AF spin structure is modulated by some func- 
tion f(x, y), the hyperfine field at position {x, y) is given 
by 

H{x,y) = ^,o{A,fix,y)-B[f{x + a,y) (3) 
+f{x -a,y) + f{x, y + a) + f{x, y - a)]} , 

whereas = 38 kOe/Ms and B = 42 kOe/nB are the 
hyperfine coupling coefficients and = 0.3/^5 is the 
maximum value of staggered polarization. Substituting 
/(x, y) given by Eq.Q, one can obtain explicitly that 
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where 



Hrn = /io [A - iBcos^inS)] 



(4) 



(5) 



It is assumed that the hyperfine magnetic fields are 
parallel to the Cu02 planes^, and that the quadrupolar 
part of the Hamiltonian is axially symmetric with respect 
to the crystal axis perpendicular to the Cu02 planes^^ 
(c-axis). With these assumptions, the full Hamiltonian 
for one ^^Cu nuclear spin has form 



n = 



2 



InKR ■ I 
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where h is the Plank constant, I is the spin 3/2 opera- 
tor, H is the hyperfine field, 7„ is the giromagnetic ratio 
of 63Cu nuclei (equal to 11.285 MHz/Tesla ), and vq 
is the quadrupolar parameter having average values^ 
34.5 MHz and 41.8 MHz for for the A- and B-hnes of 
the ^^Cu NQR spectra. The relative spectral weight of 
the B-line with respect to the A-line was extracted from 
the high-temperature NQR spectra of Hunt et al^ to be 
0.07. In the calculation, the Gaussian distribution of vq 
around the above mentioned average values was assumed 
with root-mean-squared deviation 2.8 MHz25^ 

For a given value of quadrupolar parameter i'q and 
the hyperfine field H, Hamiltonian ^ generates four 
levels and six transition frequencies between those lev- 
els (shown in Fig.[5|). 

The contribution of each transition to the NQR spec- 
trum was weighed by the square of the transition fre- 
quency, and then by the isotropic average of the transi- 
tion probability 

P„,„, = |(n'|I.Hi|n)|2, (7) 
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FIG. 6: [Color online] Distributions of hyperfine fields and 
NQR spectra corresponding to (a) commensurate diagonal 
grid superstructure shown Fig. [T] (b) incommensurate version 
of the same grid superstructure, and (c) disorder-broadened 
modification of case (b). Solid lines represent theoretical cal- 
culations. Circles represent 350 mK experimental NQR spec- 
tra from Ref.— . The relative scale of experimental and the- 
oretical NQR spectra is adjusted to give the best agreement 
between them. The inset in figure (c) shows the Fourier trans- 
form of the hyperfine field distribution from the same figure 
thus giving an idea of the amplitude of the ^SR asymmetry 
oscillations. 



where Hi is the radio frequency field (^^100 G)(See 
though the discussion of this factor in the next section.) 

When the grid structure is commensurate, it gener- 
ates only a finite discrete set of different hyperfine fields 
{Hi}. This discrete distribution was broadened to a set 

of Gaussians exp— '^^^^i^ with Ui — 0.0477j^^. 

The commensurate diagonal grid modulation was im- 
portant to test, because it implies a larger number of 
inequivalent nuclear sites, in comparison to the commen- 
surate grid considered by Hunt et al. However, as shown 
in Fig. (nja), theoretical NQR spectrum corresponding to 
the diagonal commensurate grid modulation still has too 
much structure in disagreement with experiment. 

In the case of incommensurate periodic grid modula- 
tion not considered by Hunt et al, the distribution of 
local fields was obtained in the following way: 

Grid superstructure is incommensurate, when the pe- 
riod of modulation along each of the two lattice directions 



is an irrational number [in the units of lattice period] . In 
this case, the two modulation phases are ergodic on the 
lattice, which means that, when taken modulo 27r, both 
of them uniformly fill interval (0, 27r). This implies that 
the distribution of the hyperfine fields can be computed 
from the modulation function ^ by assuming that spins 
uniformly fill the space instead of being attached to a 
discrete lattice. 

The resulting distribution function can be found ana- 
lytically. For H < Hm, it has form 



P{H) 



Arccos ( 1^ 



H„ 



(8) 



H 



where function F[£^,ri) is the elliptic integral of the first 
kind: 



F(^,7?)= / ^= 

Jo \/l — rii- 



(9) 



yysm u 



Here ^ and rj are some function arguments, and u is an 
integration variable. For H > Hm, P{H) = 0. 

In order to generate a specific distribution P{H), S ~ ^ 
was substituted into Eq.([5]), which gives H„i = 3.16 T. 
Strictly speaking S — ^ corresponds to a commensurate 
modulation but the value of the small irrational number, 
which should be subtracted from | makes very little dif- 
ference, once the incommensurate modulation sets in and 
formula ([8|) becomes applicable. 

The shape of P{H) is shown in Fig. [G^b). It has a 
Van Hove singularity at if = 0, which is associated with 
the fact that H{x,y) given by Eq.H]) has a saddle point 
at a; = 0, y = 0. This singularity, however, does not 
lead to a pronounced feature in the NQR spectrum, be- 
cause, at small values of H, the transition frequencies are 
dominated by the value of i'q, and, therefore, the NQR 
spectrum is shaped by the non-singular distribution of 
i/Q. In contrast, the step feature of P{H) at H = Hm 
falls in the range, where the value of H begins control- 
ling the values of transition frequencies. This leads to 
the [rounded] steps in the NQR spectrum at frequencies 
corresponding to H = Hm (see Fig.[S]). 

Overall, as shown in Fig. [H^b), the incommensurate 
grid modulation gives a better agreement with experi- 
ment in comparison with the commensurate one, but still 
with features in the form of steps absent in experiment. 
(A closer examination of experimental spectrum, in fact, 
reveals the hints of steps in the positions expected theo- 
retically.) 

In cuprates, there is no compelling reason to expect a 
perfectly uniform incommensurate spin density wave (ID 
or 2D). Instead, one should rather expect incommensu- 
rate periodicity only on average with local pattern dis- 
torted either by disorder or by the McMillan type of dis- 
commensurations in the commensurate superstructure^. 
A straightforward way to incorporate both effects is to 
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broaden the hyperfine field distribution obtained for tfie 
perfect incommensurate modulation. 

As shown in Fig. [6l^c), only a minor degree of broad- 
ening — 10 per cent of the original distribution width 
— is sufficient to eliminate the essential disagreement 
between the theory and the experiment. The remaining 
disagreement should be well within the uncertainty of the 
theoretical scheme and the experiment (see the discussed 
in the next section). 

The distribution of hyperfine fields in Fig. [D^c) was ob- 
tained by simple convolution of the distribution shown in 
Fig. b) with a Gaussian exp—-^, where ad — O.lHm- 
This uniform broadening was chosen over the one propor- 
tional to H [as done in Ref.^^ and for Fig. [Ha)] just for 
the sake of simplicity. In general, one should expect to 
have both 77-independent and i/-dependent terms. How- 
ever, for the present purpose, the number of significance 
is (Jd{Hm), while other differences in the choice oiad{H) 
lead only to minute differences in the NQR spectra. 

Ten per cent broadening is a very reasonable value the- 
oretically. It is expected, for example, in the case, where 
a grid line fluctuates between two adjacent diagonals on 
the lattice, which leads to the fluctuations in the size 
of the grid supercell by a/y/2. Given the approximate 
size of the supercell 4a\/2, the above fluctuation amounts 
to 12.5 per cent fluctuation in the size of the supercell, 
which, in turn, can easily lead to a 10 per cent fluctuation 
in the value of Hm within the supercell. 

Importantly, the amplitude of spin modulation used for 
the present calculation (/ig — 0.3/is ) is consistent with 
the frequency of the /iSR asymmetry oscillations-^^. One 
can further observe that the shape of the distribution of 
hyperfine fields seen by nuclear spins should not be too 
different from the shape of the distribution of magnetic 
fields experienced by muons. (Simulations in Ref.^^ in- 
deed suggest so away from the close vicinity to the zero 
value of the field.) Therefore, by Fourier transforming 
the distribution of the hyperfine fields, one can obtain 
a rough estimate of the amplitude of the /iSR asymme- 
try oscillations. Such a Fourier transform is shown in the 
inset of the right panel in Fig. [HKc) . It exhibits clear oscil- 
lations with the amplitude roughly consistent with but 
smaller than the /xSR observations. (The frequency of 
/iSR asymmetry oscillations should, of course be rescaled 
with the strength of muon-electron magnetic coupling.) 
It should be mentioned though, that, in the absense of 
NQR data, the the smallness of the oscillation amplitude 
was considered by Kojima et al^ as an argument against 
the grid interpretation. 

To summarize, none of the propositions considered by 
Hunt et alM- came close to the kind of simultaneous 
agreement with the NQR and /xSR results obtained in 
this section for the disordered incommensurate grid mod- 
ulation. 



V. DISCUSSION 

The significance of the findings presented in the pre- 
vious section should be viewed in the context of uncer- 
tainties of the theoretical framework, within which these 
findings were obtained. 

An important quantitative check of that framework is 
the fact, that it correctly predicts the experimental posi- 
tion of the 31 MHz peak. It is the check of the correctness 
of Hamiltonian ©. More specifically, since the value of 
vq is known from the high-temperature measurements, 
the above agreement amounts to a combined consistency 
check of the assumptions that the Hamiltonian is axially 
symmetric and that local hyperfine fields are parallel to 
the Cu02 planes. 

One uncertainty of the treatment in Section IIVI is as- 
sociated with the assumption of completely random dis- 
tribution of A- and B-type nuclei in the inhomogeneous 
background. One can object, that charged grid lines are 
attracted to dopant atoms, and hence the nuclei of B- 
type (attributed to the vicinity of dopant atoms) see quite 
different environment than A- type nuclei''^. It may hap- 
pen, for example, that small amplitude fluctuations of 
grid lines wipe out the signal from B-nuclei, while leav- 
ing a signiflcant signal coming from A-nuclei. It turns 
out, however, that, if only A-type (majority) nuclei are 
included in the calculation for the disorder broadened in- 
commensurate grid then the agreement between the the- 
ory and the experiment becomes even better than the 
one shown in Fig. [Sic). 

A related difficulty is that the charge modulation, as it 
becomes stabilized, creates extra electric field gradient, 
which can lead to an extra broadening of the distribu- 
tion oivQ. A priori, the relative value of this broadening 
in comparison with the one due to the intrinsic factors 
(such as charged dopants or crystal imperfections) is dif- 
ficult to estimate. However, it is also difficult to see, 
why the modulation-induced broadening would be much 
larger than the "intrinsic" one, known from the width 
of the high temperature NQR spectra. If it is not much 
larger, then the related uncertainty is not very signifi- 
cant. 

Another important element missing from the descrip- 
tion in Section IIVI is what happens in the system as a 
result of the 7r/2 and tt radio- frequency (rf) pulses of the 
spin echo sequence. The computational procedure of Sec- 
tion [IV] assumes that the 7r/2 and tt pulse conditions are 
perfectly satisfied for every transition of every spin. In 
reality, this cannot be done simultaneously for all spins 
and all transitons, because (i) in powder samples (used 
in Refi^), the polarization of the rf field has a random 
orientation with respect to the lattice directions; and (ii) 
even, if it were a single crystal, the rf field at a given 
frequency affects simultaneously several NQR transitions 
and couples to those transitions through different values 
of matrix elements (which typically vary by a factor of 2). 
The difference in the values of matrix elements implies 
that the pulse conditions cannot be fulfilled simultane- 
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ously for every transition. 

Even, if one were able to manipulate each transition in- 
dividually, the use of weighing factor ([7]) would be doubt- 
ful. Such a factor is standard and appropriate, when the 
resonance is detected by the steady state absorption tech- 
nique. In the pulsed experiments, however, one power of 
the matrix element |(ri'|I • Hi|n)| is absorbed into the 
duration of the rf pulses, which means that, if the 7r/2 
and TT conditions are perfectly satisfied, one should ex- 
pect a weighing factor equal to the above matrix element 
and not to its square as Eq.(I7]) suggests. When the pulse 
conditions are not perfectly satisfied, higher powers of 
that matrix element may contribute to the weighing fac- 
tor, but even then, the use of the factor ([7|) would likely 
amount to a drastic oversimplification. Expression ([7|) 
was used in Section HVl onlv to make the results directly 
comparable to the calculations of Hunt et al. 

Yet another uncertainty of this study is the lack of 
theory describing the wipeout of the NQR signals on the 
timescale of rf pulses. In experiment of Rcf.^-, the dura- 
tion of the rf pulses was chosen to be the one that maxi- 
mizes the spin echo response. The wipeout effect would 
suppress the intensity of the echo, if the tt/2 and tt condi- 
tions required too long pulses, and, therefore, could factor 
into the experimental definition of the pulse conditions. 

All unknown factors mentioned above can add up to 
a quantitatively significant weighing factor on the top 
of the theoretical calculation described in Section. IIVI 
This extra factor, however, is not expected to have strong 
frequency dependence and, therefore, would not mask 
the presence of the multiple peaks associated either with 
discrete or with singular distributions arising from ID 
modulations of hypcrfine fields. 

It is, therefore, the absence of the extra peaks beyond 
the one at 31 MHz that is the best indication against ID 
stripe scenarios. 

The non-observation of multiple NQR peaks also con- 
stitutes an indication against spin spirals^^ in 1 /8-doped 
lanthanum cuprates. If the Cu spins are indeed parallel 
to the Cu02 plane, and if the NQR tensor is axially sym- 
metric with respect to the c-axis, then the spiral arrange- 
ment of spins implies that all nuclear sites are equivalent, 
which should lead to sharp multiple NQR lines. If any of 
the above two assumptions is incorrect then NQR spec- 



trum should still exhibit Van Hove singularities similar 
to those shown in Figdfd). 



VI. CONCLUSIONS 

The main message of this paper is that the absence of 
multiple peaks in the NQR spectra of Hunt et al. can 
be easily explained on the basis of slightly disordered 
incommensurate grid modulation. The degree of disor- 
der required for this interpretation entails approximately 
10 per cent broadening of the non-disordered NQR line, 
which would be consistent with the fact that the signa- 
tures of spin modulations are observable by /iSR. At the 
same time, all alternative interpretations of the above 
spectra in terms of ID spin modulations encounter seri- 
ous difficulties. It is, therefore, concluded that the com- 
bined //SR/NQR data favor the the 2D grid modulation 
over the ID stripe or spiral modulation static or dynamic. 

Note added: After the present manuscript was sub- 
mitted, this author became aware of the preprint 
by Christensen et ali^, reporting significant new ex- 
perimental results of magnetic neutron scattering in 
Lai.48Ndo.4Sro.i2Cu04. These results were interpreted 
by the above authors as strongly supporting to the ID 
stipe interpretation. The present author has expressed 
his reservations on this matter in Ref4^. 
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APPENDIX 

Figure[3]shows relative peak intensities for charge mod- 
ulation ^ . These intensities were obtained by the direct 
expansion of modulation ^ [with a — I] into Fourier 
harmonics: 
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The terms grouped in parentheses (...) give the same peak positions in the first Brillouin zone. The relative peak 
intensities exhibited in Fig. [3] are the squares of the Fourier amphtudes appearing in the above expansion. 

The foUowing examples illustrate that the leading character of the peak at (±^,0) and (0, depends only on 
the presence of the sin^ factors, which suppress the charge density near the grid nodes. 
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